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Abstract 

Consider a surface, enclosing a fixed volume, described by a free-energy depending 
only on the local geometry; for example, the Canham-Helfrich energy quadratic in the 
mean curvature describes a fluid membrane. The stress at any point on the surface 
is determined completely by geometry. In equilibrium, its divergence is proportional to 
the Laplace pressure, normal to the surface, maintaining the constraint on the volume. 
It is shown that this source itself can be expressed as the divergence of a position- 
dependent surface stress. As a consequence, the equilibrium can be described in terms of 
a conserved effective surface stress. Various non-trivial geometrical consequences of this 
identification are explored. In a cylindrical geometry, the cross-section can be viewed as 
a closed planar Euler elastic curve. With respect to an appropriate centre the effective 
stress itself vanishes; this provides a remarkably simple relationship between the curvature 
and the position along the loop. In two or higher dimensions, it is shown that the only 
geometry consistent with the vanishing of the effective stress is spherical. It is argued 
that the appropriate generalization of the loop result will involve null stresses. 
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1 Introduction 

The energy of a soap bubble is proportional to its surface area. The bubble is round because 
this is the shape of the surface, enclosing a fixed volume, which minimizes area [T]. The shape 
of a closed surface described by an energy which depends in some more general way on the 
geometry need not be spherical. A striking example is provided by biology; the fluid membranes 
in our cells are described with uncanny accuracy, on optical scales, by the Canham-Helfrich 
(or Willmore) energy quadratic in the mean curvature of the membrane surface E]- ' n 
addition to the volume the area is also fixed. There may also be other constraints. Under these 
circumstances, determining the equilibrium shape, never mind addressing questions of stability, 
becomes a sutble problem. We know, however, thanks to a series of landmark numerical studies 



1 



in the nineties that very different equilibrium shapes may be consistent with given values of these 
constraints and that large changes in shape may be brought about by tweaking the constraints 
slightly. This work is nicely summarized in references [4 5J and |6J. Representative examples of 
recent work are provided by [7J and [8J. Our current theoretical framework, by constrast, remains 
woefully undeveloped. Even if we artificially limit ourselves to axially symmetric configurations, 
where it is possible to boast a level of analytical control, there is a significant discrepancy 
between the wealth of numerical data that has accumulated and our current ability to interpret 
the patterns within it — even qualitatively — on the basis of the underlying model. Some new 
element needs to be introduced into our theoretical framework if this imbalance is to be set 
right. 

It is well known that is it possible to describe the equilibrium shape of the surface in 
terms of the conservation of a stress tensor. This problem was first approached, a la Gibbs, 
thermodynamically |9][TU|; such an approach, however, never really exploits the fact that the 
energy depends only on geometrical degrees of freedom; when the scaffolding is removed, what 
is left is geometry. The appropriate setting has been described in [TT][T2]. Recently, Lomholt 
and Miao have shown how the two approaches can be reconciled [T3]- From the geometrical 
point of view, the existence of the conservation law is a consequence of the fact that the free 
energy is invariant under translations. Noether's theorem identifies the stress as the conserved 
current associated with this invariance. Furthermore, the invariance of the geometrical free- 
energy under reparametrizations of the surface implies that the stress is determined completely 
by the surface curvature and its derivatives. Knowledge of the local membrane geometry will 
always be sufficient to reconstruct the stress underpinning it. This contrasts with the behavior 
in an elastic solid where the stress is proportional to first derivatives and the geometry itself is 
not always enough to determine the state of stress [14J. 

Internal constraints on the geometry will add stress to the membrane. In the case of a fluid 
membrane, as we have mentioned, the area is fixed; the two layers of a bilayer may also have 
different areas. The latter asymmetry is captured by a constraint or a penalty on the total mean 
extrinsic curvature of the surface [15J. These constraints are accommodated using Lagrange 
multipliers. 

Whereas surface constraints add to the stress due to bending energy, a constraint on the 
enclosed volume acts as an external force on the membrane. This is the well known Laplace 
pressure counteracting the tension in a soap bubble which acts normally to the surface. In this 
note, it will be shown how to incorporate this external pressure into a divergence-free effective 
surface stress. This makes it possible to treat the volume constraint on an equal footing with its 
surface counterparts. Geometrically, this is because the surface normal along which the Laplace 
force acts can always be expressed as a surface divergence: a fact which is itself a direct 
consequence of the translation invariance of volumes. Curiously, the translational invariance of 
the stress gets broken in the construction; the effective stress depends explicitly on the position. 

The interesting point is that this unprepossessing mathematical identity permits us to look 
at the equilibrium in a new light. We first examine a model describing cylindrical geometries in 
terms of this effective stress tensor. The cross-section of a cylinder perpendicular to its axis is 
a closed planar loop; the free-energy on this geometry describes one of Euler's elastica subject 



2 



to a constraint on the enclosed area on the plane. This system was studied in some detail in 
[16J where it was dubbed, with a nod to Euler, Elastica Hypoarealis. While it is well known 
that the Euler-Lagrange equation can be integrated to provide a quadrature for the curvature, 
this system exhibits an additional level of integrability which is not evident in its Euler-Lagrange 
equation: the curvature at any point detemines completely the geometry of the loop at that 
point. Such an identity is surprising because one would expect to require two integrations of 
the curvature to reconstruct the corresponding loop whereas none is needed. In this paper, 
it will be shown that this is no accident. In fact, it is because the conservation law can be 
integrated: modulo a translation centering the loop at the origin, the effective stress vanishes 
everywhere on the loop. The identity between curvature and position pointed out in [16J follows 
as an immediate consequence. Given the curvature, it is possible to trace the loop. One has 
not, of course, integrated the shape equation: one still needs to integrate the quadrature for 
the curvature as a function of arc-length in terms of elliptic integrals. But one does possess 
a relationship between curvature and the configuration which has sidestepped completely the 
need to evaluate these integrals. 

The construction of the identity for a loop is loosely analogous to the Wulff construction 
of the equilibrium shape of a crystal from a knowledge of its anisotropic surface tension [T7j . 
Is there an analogue for genuine two-dimensional membranes? Just as the equilibrium loop 
possesses a centre of 'gravity', there also will be one for a membrane. The relationship between 
curvature and the location with respect to this centre originates in the explicit spatial dependence 
of the effective stress. In a spherical equilibrium, it is simple to check that this stress vanishes 
when the origin is translated to the centre of the sphere. This is, in itself, a non-trivial identity 
concerning the balance of internal forces in the spherical equilibrium. It is, however, also simple 
to prove in the case of a fluid membrane that spheres are the only compact geometries possessing 
this property. 1 One does not appear to have new to say about the shape of membranes. 

There is, of course, a good reason why a straightforward analogue of the loop identities does 
not exist in higher dimensions. On the loop, a divergence-free (one-dimensional) effective stress 
is trivially constant, and this constant can be set to zero by centering the loop appropriately. In 
two-dimensions, however, there is enormously more freedom. The analogue of a constant vector 
in the plane is now a null stress, divergence- free by construction. The null stresses on a surface 
are in one-to-one correspondence with spatial vector fields. Furthermore, consistency with the 
symmetries of the fluid membrane stress tensor will place a constraint on this vector field; one 
finds that the admissible vector fields are the generators of area-preserving deformations of the 
surface. Such stresses cannot generally be eliminated by a simple translation. As a consequence, 
the higher dimensional generalization of the result for loops will involve null stresses. It is too 
early to know if this going to be tractable. However, given the enormous stakes involved (and 
not just for fluid membrane physics) this is certainly worth exploring. 

1 1 n the case of the Clifford torus, one can check explicitly that the effective stress is non-vanishing. 
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2 First order variations, Euler-Lagrange derivatives, translation in- 
variance and the stress 



The membrane will be described as a parametrized surface, £ a — > X(£ a ). The X = (X 1 , X 2 , X 3 ) 
are three functions of two parameters £ a , a = 1,2. The metric tensor induced on this surface is 
given by g ab = e a -e b ; the corresponding extrinsic curvature is K ab = e a -d b n, where e a = <9 a X, 
a = 1,2 are the two tangent vectors to the surface and n is the unit normal. V a represents 
the covariant derivative compatible with g ab . The Gauss-Weingarten equations, describing how 
these vectors change as one moves across the surface, are V a e& = — K ab n and V a n = K a b e b . 
Indices are raised and lowered with the metric. When one looks at cylinders, the surface of 
interest will be a planar curve. Higher dimensional generalizations are also of potential interest; 
because it is straightforward to frame the discussion in terms of D-dimensional surfaces this 
is what will be done unless explicitly indicated otherwise. Geometrical background material is 
provided in jTHJ [T2|. A nice tutotial tailored to membrane geometry is provided by reference 

Consider any geometrical integral of the form 

H[X]= JdAH(g ab ,K ab ), (1) 

defined on a patch of surface £ with a boundary <9£. H will represent the free energy of the 
surface. H is a function of the metric and the extrinsic curvature. Thus H will be a functional 
of X and it will be invariant under reparametrizations and Euclidean motions of the ambient 
space. The area element induced on the surface is dA = ^det g a bd 2 £. For simplicity, gradients 
of K ab are not considered; they have, however, been considered in some detail elsewhere [21J. 

Recall how the stress tensor corresponding to H gets constructed. This involves examining 
the response of the free energy to a small deformation of the surface. This deformation is 
described by the infinitesimal change in the embedding functions: X — > X + 5X. The response 
of H decomposes naturally into two parts: a bulk response as well as a boundary one; the latter 
is a surface divergence originating in derivatives of the deformation in the bulk which have been 
consigned to the boundary. One finds that 



5H = J dASn- 5X- J tL4V a [F • <JX] 



dAVn 



H ab -V b 5X . (2) 



Here £n is the Euler-Lagrange derivative of H with respect to X; because of reparametrization 
invariance of H, only the projection of <5X onto the normal contributes to the bulk variation. 
f a is the surface stress tensor associated with if and H ab = &HjdK ab . It has been shown 
elsewhere [TT] that 

£ = - V a V b H ab + K ac K c b H ab + UK ■ (3) 

and 

f a = (H ac K c b - Hg ab )e b - V b H ab n . (4) 
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When H is translationally invariant it is possible to rewrite its Euler-Lagrange derivative as a 
surface divergence. The argument of this divergence is the stress. This follows from Eq.fl^l): 
consider the effect on the surface of a translation, <5X = a, where a is a constant vector. 
The derivative term appearing on the second line in Eq.(J2J vanishes; translational invariance, 
SH = 0, then implies that 

a ■ J dA [8 n - V„f a ] = . (5) 
Because the region of integration itself is arbitrary, the local statement 

8n - VJ a = (6) 

follows; 8 n is a surface divergence. This is a useful identity to remember. 
In particular, a fluid membrane is described by the Canham-Helfrich energy 

H Q = l ~jdAK\ (7) 

where K = g ab K a i,. It's mathematical properties, notably its conformal invariance, were exam- 
ined long beforehand [3J. Suppose that both the area A and the integrated mean curvature 

M = J dAK (8) 

are fixed. The free-energy describing the surface is then given by 

H = H + (3(M- M ) + a(A - A ) , (9) 

where [3 and a are two Lagrange multipliers enforcing the constraints. The constants M ,A 
are the (for our purposes) irrelevant fixed valued of M and A. 



3 Equilibrium with constrained volume 

Consider now the problem of determining the equilibrium of a closed membrane subject to a 
constraint on the enclosed volume. Introduce a new energy functional H related to the surface 
free energy H by 

H = H-P(V-V ). (10) 

The constant Laplace pressure P is a Lagrange multiplier forcing the volume to assume some 
fixed value, Vo- In equilibrium, the Euler-Lagrange equation 8 = P is satisfied (8 is the 
Euler-Lagrange derivative of H). Using Eq.fJSJ), this equation can be cast as a conservation law: 

V/ = Pn. (11) 

The external pressure provides a source for the stress in the membrane. It is always possible, 
however, as we will now demonstrate to incorporate this particular source into an effective 
conserved surface stress. 
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4 Cone volume and the effective stress 



The stress was identified in section 2 by examining the boundary behaviour of the energy under 
deformation. One way to place the volume V appearing in the constrained energy given by 
Eq. (ITU|) on an equal footing with the surface invariants appearing in H is to associate a notion 
of volume with patches of the surface. But this is easy to do: fix an origin, and construct 
the cone C with a base which consists of the patch and its apex positioned at the origin. By 
examining how the volume of this cone behaves under deformations of the surface patch, we 
will see that it is possible to associate an effective surface stress tensor with the volume term 
appearing in Eq. ()10|) . It it then straightforward to recast the equilibrium condition under the 
external Laplace pressure in terms of a divergence-free effective surface stress. 
First consider a closed surface. The volume of space it encloses is given by 

V = — *— [ dAxi- X. (12) 
D+l J K ' 

This is a direct consequence of the elementary Euclidean space identity for the position vector 
x, div ■ x = D + 1 and the use of the divergence theorem. While it is possible to express V as 
a surface integral, it is not possible to express it in the form (pfj) . Despite this, as we will now 
show, an analogue of Eq.fl^J) does hold. 

Suppose that the integral on the rhs of Eq. (IT2"l) is restricted to some surface patch. It is 
easy to see that the integral then represents the volume of the cone C that was described at 
the beginning of this section (see |19J). 

Consider now the response of the volume of this cone to a deformation of the surface 
X — > X + <5X. This will point us towards an extremely useful identity for the normal vector. 
One finds that 

5V = J dAn- 5X- J dAV a (tf ■ SX) , (13) 

where 

f a = (e a • X) n - (n ■ X) e a = X x (n x e a ) . (14) 

The second expression appearing on the rhs of Eq. ()14|) is valid for two-dimensional surfaces. The 
derivation is straightforward; the details are provided for the interested reader in an appendix. 
In particular, the well known result that the Euler-Lagrange derivative of volume with respect 
to X is n is reproduced. On a closed surface, the second term in Eq. (ll3|) vanishes so that the 
familiar area times normal displacement formula for the shell volume SV follows. 

The stress-like object fg we have introduced has a very nice property: its surface divergence 
is proportional to n: 

n=^V a f Q . (15) 

Curiously, the fact that n is a surface divergence appears to have gone unnoticed. The value 
of Eq. (ll5|) is that it provides a relationship analogous to Eq.flf)]) for the volume with 8 = 1: 
ifi/D is identified as the effective surface stress associated with the volume term in the energy. 
The derivation of Eq. (IT5|) is elementary: note that the curvature terms, which originate in the 
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Gauss- Weingarten equations for the surface, cancel when the two terms are summed. This is 
consistent with the fact that the normal depends only on first derivatives of X. 

A consequence of this identification is that the Laplace force on a patch of surface can be 
cast as a closed line integral along its boundary: 

pjdAn = -^s(Xxt). (16) 

Here t is tangent to the boundary curve and s is arc-length. 2 In this form, it is clear that the 
Laplace force on the patch depends only on the boundary geometry, a fact that is not at all 
obvious when cast as an integral over the patch. 

Note that the effective stress is not f"o/(-D + 1) as a casual comparison of Eq.(lT3|) with 
Eq.(0) might suggest; this is because the cone volume, unlike the volume of a closed surface, is 
not translationally invariant if the apex of the cone is fixed. Indeed under a constant translation 
a, X — > X + a, and it is clear from Eq. ()12|) that 

SV = -z— — a • / dAn, (17) 
D + 1 J v ' 

a result consistent with Eq.(lT3|) with the replacement of 5X by a. 

The structure of the boundary term in Eq.(ll3p is of interest. On a two-dimensional surface, 
it is possible to cast this equation in the form 



5V 



J dAn-6X + -fds(txX)-5X. {It 



In this form, it is clear that the boundary term depends only on the boundary curve and its 
deformation; in particular, it is independent of what the patch gets up to in its interior. 



5 Laplace pressure as equivalent surface stress 

The identification (fT!)j) makes it possible to recast the equilibrium condition under an external 
Laplace pressure in terms of a divergence-free effective surface stress. Introduce the effective 
stress tensor, 

fa = f a_ P f a (lg) 

where fg is given by Eq. (IT4|) . Using Eq. (IT5|) . it is possible to recast the equilibrium condition 
Eq. flTTj) in the form 

v a f° = o. (20) 

i a unlike f a is divergence- free in equilibrium. In an isolated membrane, it follows from Eq. ()20|) 
that 

<fdslj a = (21) 

2 We have used the fact that for right-handed triad {1, t, n}, t = n x 1, where 1 = l a e a is its outward pointing 
normal to the curve lying on the surface. Any two orthonormal vectors specify a third; in particular, the two 
normals n and 1 to the curve fix its tangent. 
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along any contractible curve which lies on the surface. With the original stress tensor, the 
corresponding integral depends on the behaviour of the normal vector throughout the enclosed 
patch. Using the effective stress tensor, the volume constraint is treated like the other surface 
constraints. 

If an external force acts on the patch of membrane enclosed by the curve, the line integral 
appearing in Eq.(l21|) will determine this force. In this context, it is often possible to deform the 
contour appropriately to take advantage of any symmetry the configuration might possess. This 
approach was developed in [22] to examine surface mediated interactions. Its generalization to 
incorporate a pressure difference across the membrane is also possible. 

The effective stress f a has one very striking feature we have yet to comment upon. Unlike 
f a , f a is not translationally invariant. Its origin is the explicit X dependence appearing in fg 
given by Eq.(ll4|). Under a translation X — > X + a, one finds that fg — > fg + hg, where 

hj = ax(nxe"). (22) 

This shift should leave no physical trace. In the next section, this issue will be examined in the 
appropriate setting. 



6 Null stresses 

It is clear that there is an inherent ambiguity in the definition of any divergence-free effective 
stress f°; this is because one is always free to add to f a any stress of the form 

h a = V b A ab , (23) 

involving an antisymmetric potential A ab = — A 6a . (In the case of a curve, this ambiguity 
degenerates into the freedom to add a constant vector.) It is straightforward to show that for 
h a constructed according to Eq . ()23|) . one has V a h a = or h a is conserved. We will refer to 
such an h" as a null stress. In the language of differential forms h a is exact |23J. 

On a two-dimensional surface, the anti-symmetric potential A ab factorizes: it must be 
a product of the two-dimensional Levi-Civita anti-symmetric tensor e ab and a space vector 
potential A: A ab = e ab A. Because e ab is covariantly conserved (V a Q, c = 0), it is possible to 
express the corresponding null stress in the simple form [TTJ 

h a = e ab _ ( 24 ) 

How various ambiguities in the definition of f a get captured in h a are discussed very nicely in 

Suppose now that we add an h a constructed this way to the stress; its contribution to 
the force on any patch is given by the line integral of l a h a along the closed boundary curve. 
Because h a is exact, this integral vanishes. This is well known. To see it explicitly when D = 2, 
express e ab in terms of the components of the unit tangent t = t a e a and the normal 1 = l a e a , 

e ab = lat b _ t b t a SQ that 

j> ds l a h a = - ids t a V a A . (25) 
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But t a V a = d/ds so that the closed line integral vanishes. It would appear that the notion of 
a null stress is not very interesting; as we will see this is not the case. 

We are now in a position to comment on the behaviour of under translations. It was 
seen that under a translation X — > X + a, fg transforms by fjj — > + hg, where hg is defined 
by Eq. ()22|) . It is easy to see that is, in fact, null. Indeed, when D = 2, 

h a = e ab a xe b = e ab V b [si x X] , (26) 

which is of the form ()24|) with the identification, A = a x X. 

While a change of origin has no physical significance, the value of f a at a point will depend 
on the choice made. It may be useful to centre the origin to reflect the symmetry of the 
configuration. For example, in a sphere it would be perverse to place the origin anywhere but 
at the centre of the sphere unless there is a convincing reason to do otherwise. In the following 
section we will see that the explicit X dependence which appears in f a has a very important 
role. 



7 An identity for rigid loops enclosing a fixed area 

A strong hint that casting the conservation law in terms of an effective stress tensor might be 
useful is provided by examining a membrane which forms an infinite cylinder. The cross-section 
of the cylinder perpendicular to its axis is a closed planar loop. The free energy for the cylinder 
ascribes rigidity to the loop. Its length L is fixed and it encloses a fixed area A. Consider now 
a curve embedded on the plane parametrized by arclength s, X = X(s); a prime will denote 
a derivative with respect to s, so that t = X' is the unit tangent vector. The curvature k is 
defined by t' = — kn and n' = kt, where n is the unit normal. The energy penalty associated 
with bending is given by 

H = J dsk 2 . (27) 

To facilitate comparison, we will use the same normalization for the bending energy introduced 
in [16J (without the factor of a half), but adopt a more conventional soft-matter notation for 
the multipliers. Now the energy functional capturing the constraints on L and A is 

H = H + aL- PA. (28) 

Consider the problem of determining the equilibria of this loop. The Euler-Lagrange equation 
is 

2k" + k 3 -ak + P = 0, (29) 

and it has the first integral 

h' 2 + V(k) = E, (30) 

where 

V(k) = -k A - —k 2 + —k , (31) 
v i g 4 2 > v ) 
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and E is a constant. The closure of the loop will quantize E. The quadrature determines 
the curvature k implicitly as a function of s in terms of elliptic functions. One would expect 
to require two integrations of the curvature to reconstruct the corresponding loop; remarkably, 
none is needed. This was shown in [T§]. Here, it will be shown how this comes about using the 
conservation of the effective stress. 
The tension in the loop is given by 

T = (k 2 - a)t - 2k'n , (32) 

and it satisfies 

T' = Pn. (33) 

For a planar curve (D = 1), Eq. f)15j) simplifies n = [(X ■ t) n — (X • n) t]'. The equilibrium 
condition Eq.(l33p may thus be recast as an identity for a space vector 

[X ■ n + p-\k 2 - a)] t - [X • t + 2P- 1 k'} n = C , (34) 

where C is some constant vector along the loop. This constant is the analogue of h a discussed 
earlier in the two-dimensional context. However, it is always possible to translate the origin so 
that X • t = when k! = 0, which gives C = 0. Eq. (l34|) now provides two identities, one for 
each projection 

Xn = p-\a-k 2 ) 

X t = -2P- l k' . (35) 

The V identities which were identified in [TH] are reproduced in a very direct way. The latter 
of the two identities implies that 

X 2 = X 2 - 4P -1 Jfe , (36) 

where X is a constant. Completeness of the basis vectors t and n together with Eq . ()30|) can 
now be used to express Xq directly in terms of the constant E appearing in the quadrature: 

X 2 = P~ 2 (8E + a 2 ). (37) 

Once the spectrum of E has been determined, the curvature at any point detemines completely 
the loop geometry at that point. No integrations are necessary. While it is well known that 
the Euler-Lagrange equation can be integrated to provide a quadrature for the curvature, this 
additional level of integrability is not evident in the Euler-Lagrange equation. In a beautiful 
paper, Joel Langer mentions the existence of such an identity somewhat obliquely as a property 
of the torsion-free limit of a certain integrable model for space curves [24J. 

It should be remarked that there is considerably greater freedom associated with an infinite 
or open elastic curve due to boundary conditions. No simple relationship between curvature 
and position is known. Nor is it expected. The configuration will depend sensitively on the 
boundary conditions (see, for example, |25j). Isolation is important. 
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8 Vanishing effective stress implies spherical symmetry 

It is not unreasonable to inquire if there exist identities between curvature and position in 
solutions to V a f a = analogous to those for the planar elastic loop. We begin with a search 
for solutions satisfying the straightforward analogue to a vanishing C in Eq. (]3*4|) . P = 0. 

Recall that translational invariance and reparametrization invariance together imply a set of 
integrability conditions on the components of the stress |TT|. These conditions can be isolated 
explicitly by expanding f a with respect to the basis {e a ,n} 

f fl = A + fn. (38) 

Here f ab is a surface tensor; f a is a surface vector. The projections of Eq.(J13) perpendicular 
and parallel to the surface then read 

V a f a -K ab f ab = £, (39) 
V a f ab + K h a r = 0. (40) 

The (normal) Euler-Lagrange derivative is captured completely by the normal projection of the 
divergence; the projection process dismantles the natural surface divergence of f a . 

The remaining D identities (14Q*|) provide consistency conditions on the stress; they are 
completely independent of £; thus, they are also independent of the equilibrium. In particular, 
if H is a sum of terms, as it is for our constrained fluid membrane, Eq. ()40|) will hold for each of 
them separately. In Eq. (J39"|) . both the tangential stress f ab and the Euler-Lagrange derivative 
£ act as a source for the normal stress. In Eq. (l40|) . the normal stress provides a source for 
the tangential stress. Neither stress alone will generally be conserved. This coupling between 
normal and tangential stresses should be contrasted with the simple intrinsic conservation law 
for the stress T ab associated with a field coupling to the intrinsic surface geometry, V a T ab = 0. 

Using Eq.flU), it is straightforward to identify the stress in a fluid membrane with an energy 
given by Eq.(J9]). One finds for the surface tensor f ab and the vector f a : 

f ab = K(K ab - jg ab ) + (3{K ab - Kg ab ) - ag ab 
fa = _ v a R 

The corresponding Euler-Lagrange derivative is given by [TT1 125] 

£ = - V 2 K + \k(K 2 - 2K ab K ab ) +/3TZ + aK. (41) 

The intrinsic scalar curvature 1Z is related to the extrinsic curvature by the contracted Gauss- 
Codazzi 

K = K 2 - K ab K ab . (42) 
This setup was described in [TI] and refined in [12J. Note the following points: 

(1) The surface tensor f ab is quadratic in g ab and K ao . 
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(2) In a soap film, / 



ab 



ag ab is isotropic, with tension a; in general, however, the mechanical 



surface tension is not a. The tangential stress f ab depends not only on g ab but also on the 
local state of bending as characterized by K ab . 



(3) Because / depends locally on K ab , its eigenvectors coincide with those of K ab ; the 
corresponding eigenvalues are quadratic functions of the principal curvatures. 

(4) Umbilical points are those where the principal curvatures coincide, so that K ab = g ab K/D 
(all points on a spherical membrane are umbilical). At such a point (when D = 2), the bending 
contribution to the stress vanishes. What stresses are present are there because of competing 
constraints on the membrane geometry. 



The components of the effective stress are given by 



P 



,ab 



f ab = f ab + -(n-X)g° , 



and 



f a = f a - — (e° 



X) 



(43) 



(44) 



Various interesting scalars can be constructed using the tangential effective stress f ab given by 
(|43*|) . In the appendix, two of these scalars are shown to be surface divergences in equilibrium. 

It is simple to show that f a = in a spherical equilibrium when the origin is located at 
the centre of the sphere. To confirm this, note that on a D-dimensional sphere of radius R, 
K = D/R; in addition, on a sphere every point is umbilical, K ab = Kg ab /D, so that the 
Euler-Lagrange equation £ = P reduces to the vanishing of the following cubic polynomial in 
R: 

4 i 1 ~ d) + 2f3R + aR2 ~ pr3 = ° ' ^ 

This equation degenerates into a quadratic when D = 2 because of the scale invariance of the 
Canham-Helfrich energy. Now, if the origin is placed at the centre of the sphere, X = _Rn. It is 
straightforward to show that f ab = coincides with the identity Eq.fl45j). In addition, f a = 
because K is constant and X is orthogonal to the surface. Thus f a = 0. This is independent of 
the dimension. In particular, the scale invariance of the Canham-Helfrich energy when D = 2 
does not play a role. 

The only equilibrium configurations consistent with f a = are, in fact, these spheres. This 
follows from the vanishing of f ab . To see this, note that f ab = can always be cast in the form 



where 



K ab = Fg 



ab 



K 2 P 

— +(3K + a--(n-X) 



K + f3 



(46) 
(47) 
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In general, one can use completeness of the (orthonormal) principal vectors l l a and l 2 to provide 
a spectral expansion for K ab (Ki and K 2 are the principal curvatures) 

K ab = K x + K 2 l a 2 l h 2 . (48) 

Eq . ()46|) then implies that the principal curvatures must coincide everywhere: K\ = F = K 2 . 
What is more, they must also be constant. The simplest way to see this is to appeal to the 
contracted Codazzi-Mainardi equations, VbK ab — V a K = 0, which imply that V a F = 0. Thus 
F must be constant on the surface. The only solution is therefore spherical. Note that there 
is no corresponding constraint on the curvature of a loop (when D = 1 the Codazzi-Mainardi 
equations are vacuous). 

There are a few noteworthy alternative formulations of this non-existence result. Consider 
the following: if f a = in equilibrium, the equilibrium must be spherical and f ab = 0. 3 Note 
that if f a = 0, where f a is given by Eq.fHU), then 

K — C — — PX 2 , (49) 
2D K ' 

where C is a constant. This looks like a very direct analogue of Eq.(l36p. It turns out, however, 
that this equation is generally inconsisitent with its tangential counterpart; indeed the only 
equilibrium consistent with it is spherical. For if f a = then V a f ab = in equilibrium. But, 
on substituting the the expresssion for f ab given by Eq. (B3j) . this statement reads 

K ab [V a K-^(e a -X)] = 0. (50) 

If Eq. ()49p is now used to express K in terms of X, there follows X ■ e a = or X 2 is constant. 
Spherically symmetry is forced on us. In the next section we will explain why the naive analogue 
of the loop identities comes up short and suggest a generalization. 



9 Is there a generalization of the loop identities? 

It has been shown that a direct analogue of the loop identities does not exist unless the geometry 
is spherical. There is, however, a good reason for this failure: on the curve, the only null stress 
is constant; on a two-dimensional surface, however, any space vector field will generate a null 
stress. A more appropriate analogue is an identity of the form, f a = h a , where h a is some 
(hopefully) local null stress. Let us express h a = h ab e a + h a n. In components, 

f ab + ^(n-X)g ab = h ab , (51) 

and 

r-^(e a -X) = h a . (52) 
3 lf f a = 0, then K a bf ab = must also hold when the Eli ler- Lagrange equation is satisfied. 
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To consistently couple h ab to a symmetric f ab , h ab must also be symmetric. On one hand, this 
symmetry does not itself possess an invariant significance. This is because a null addition of the 
form t^^ will not generally have a symmetric tangential component. However, the canonical 
i a for a fluid membrane does come with this symmetry so if we are to couple h a consistently 
to f a , it had better also possess it. 

It is simple to characterize null potentials A (defined in section^ which yield symmetric h ab . 
In appendix ([0]) we show that for D = 2, h ab is symmetric when V a A ■ e a = 0. This identifies 
the vector potential A with an area preserving deformation of the surface. It is straightforward 
to construct null stresses of this form. They are also of considerable interest in their own right 
leading to unexpected connections between geometrical invariants. Work in this direction will 
be reported elsewhere. 

Whereas the symmetry of h ab is a necessary condition on h a , it is clearly not a sufficient 
one. Indeed, it is straigtforward to construct a counterexample. Suppose that h ab is a linear 
combination of g ab and K ab . The null stress h a derived from the potential A = X x n is of 
this form. It is then clear from Eq. (j51|) that K ab must be proportional to g ab . An identical 
argument to the one employed earlier kills all but the spherical solution. Unfortunately, we are 
not even close to an exhaustive list of necessary and sufficient conditions on h a . It will clearly 
be necessary to introduce some element of anisotropy into h a if spherical symmetry is to be 
broken: its symmetry will reflect the symmetry of the configuration. 

It is not, of course, obvious that Eqs. (]5ll and (|52|) are going to be tractable; they are very 
different from the Euler-Lagrange equations we are familiar with. What is clear is that it is 
worth finding out; for, if they do turn out to be tractable, the implications will be far-reaching. 

10 Discussion 

We have shown that the Laplace pressure on a closed vesicle enforcing the constraint on the 
enclosed volume can be incorporated into a conserved effective stress tensor. What might 
appear to be some mathematical sleight-of-hand does turn out to lend insight into the nature 
of the surface geometry. 

What is perhaps most intriguing about this set-up is that it suggests that, in an isolated 
vesicle, it may be possible to integrate the conservation law. In a loop, where this is straight- 
forward to do, it leads to powerful identities between the curvature and the position; control is 
provided over curvature without any need to solve the Euler-Lagrange equation itself explicitly. 
The generalization to higher-dimensional surfaces, involving the identification of the effective 
stress tensor with an appropriate null stress, is unfortunately anything but straightforward. The 
challenge is to implement it. 

It should be pointed out that our results can also be cast in the language of differential 
forms (see, for example, [27J, and in the specific context of membranes |28j). This is also likely 
to be useful. 
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A Derivation of Eq.flEU) 

To derive Eq. (lT31) . recall that the deformation SX induces the change Sg ab = e a • V b SX + e b ■ 
V a <5X in the metric, so that SdA = dAe a ■ V a <5X. The deformation in the unit normal is 
tangential, 5n = — (n • V„5X) e a . Thus 

SV = J dA [(n • X) (e a • V a 5X) - (e a • X) (n • V a 5X) + (n ■ 8X)} 

= J dA [(n ■ SX) - (f a • V a 5X)} 

= JJ~^i J dA K n ' 5X ) + ( V - f o ' 5X )] " Jj~p[ J dA V -( f o • 5X ) ' ( 53 ) 

where Eq. ()14|) has been used on the second line and all derivatives of SX have been collected 
in a divergence on the third. Now use the identity Eq.(ll5p on the second bulk term. Eq. f)13j) 
follows. 

B Identities involving e ab 

For two-dimensional surfaces, it is useful to express fg as a double cross product: 

f a = X x (n x e a ) = e ab X x e b . (54) 

Here e ab is the two-dimensional Levi-Civita antisymmetric tensor. We have e a x e b = e a fen, 
so that n x e a = e ab e b . We note that e ab e cd = g ac g bd — g ad g bc - t an d, as a consequence, 

a ac a b n ab 

fc t c y 

C Two scalars 

We examine various scalars that can be constructed using the tangential effective stress f ab 
given by (|43|). 

The trace g a bf ab is given by 

9abf ab = (l - y) K 2 + (1 - D)(3K -Da + P(n-X). (55) 
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In general, g a bf will be a divergence in equilibrium. This is a consequence of the fact, observed 
in (TTJ, that 

g ab f ab = V a (X-f a ), (56) 

when V a f a = 0. This identifies the argument of the divergence as X- f a . Upon integration of 
Eq. (|55|) over the closed surface, there follows 

(2 - D)H + (1 - D) f3M - Da A + (D + l)PV = . (57) 

When D = 2, this is the well-known virial identity associated with the scaling behaviour of 
the energy in equilibrium [5J. We also note that the vanishing of g a bf ab is consistent with this 
identity. 

In addition, K ab f ab is a divergence. This follows from the equilibrium condition, V a / a + 
K ab f ab = 0. We have 

K ab f ab = \k{K 2 - 2K) -(3K-aK+^K(n-X). (58) 

The contracted Gauss-Codazzi equation ()42"|) has been used to express the quadratic K ab K ab 
in terms of 1Z and K 2 . Using the Minkowski identity for the area [19J, 

A = — J dAK(n-X) , (59) 

one sees that the identity 

JdAK ab f ab = (60) 
is consistent with the integrated Euler-Lagrange equation with £ given by Eq. (l41|) . 



D Symmetric h and area-preserving deformations 

Note that h ab = h a e b and h ab will be symmetric if e ab h ab = 0. We now use Eq.(J2H) to express 
this condition in terms of a potential: 

e ab h ab = V a A-e a . (61) 

Thus h ab is symmetric when V a A • e a = 0. But this condition identifies the vector potential A 
with an area preserving deformation of the surface: Under the deformation <5X = A , it follows 
from the definition of g ab that Sg ab = e a ■ VfA + • V a A. The area element is proportional 
to y/g, where g = detg ab . However, 

5g = gg ab 6g ab = 2ge a -V a A. (62) 

Thus b^fg = when V a A • e a = 0. This condition can be cast in terms of components 
tangent and normal to the surface, A = A^e a + Aj_n. Then V a A ■ e a = is equivalent to 

V a A? { + KA± = 0. 
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